Anomalous angular dependence of the dynamic structure factor near Bragg 

reflections: Graphite 
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The electron energy-loss function of graphite is studied for momentum transfers q beyond the first 
Brillouin zone. We find that near Bragg refiections the spectra can change drastically for very small 
variations in q. The effect is investigated by means of first principle calculations in the random phase 
approximation and confirmed by inelastic x-ray scattering measurements of the dynamic structure 
factor S(q, cj). We demonstrate that this effect is governed by crystal local field effects and the 
stacking of graphite. It is traced back to a strong coupling between excitations at small and large 
momentum transfers. 

PACS numbers: 71.45.Gm, 78.70.Ck, 81.05.Uw, 73.21.Ac, 79.20.Uv, 82.80.Pv 



The momentum resolved and frequency dependent dy- 
namic structure factor S'(q, cj) is a fundamental quantity 
in plasma physics, nuclear physics, particle physics and 
condensed matter physics. It is important for the un- 
derstanding of many problems like e.g. electronic cor- 
relation, and it links the theory of many-body systems 
to scattering experiments like electron energy-loss spec- 
troscopy (EELS) or inelastic x-ray scattering (IXS). 

Most EELS measurements are restricted to moderate 
momentum transfer with q typically being shorter than a 
reciprocal lattice vector due to multiple scattering. IXS, 
while being less appropriate in the range of very small 
q, does not suffer from this restriction [J 0]. Modern 
synchrotron radiation sources have therefore opened the 
way to study electronic excitations at large momentum 
transfer. New phenomena can be observed in this range 
such as a periodic plasmon dispersion in magnesium di- 
boride MgB2 [2,] or plasmon-bands in Silicon [3[ . Crystal 
local field effects (LFE) 0, Q , namely the fact that an ex- 
ternal perturbation can induce a response on the length 
scale of the structure of the material, become increas- 
ingly important for large momentum transfer Q. They 
acquire particular importance for layered systems, such 
as graphite, or superlattices Electronic excitations in 
layered systems can reveal aspects that differ from the 
material in its ground state. Even for a system where 
the ground state charge density is mainly confined to the 
individual layers, excitations can induce Coulomb poten- 
tials that lead to significant interlayer interaction. One 
finds for example that plasmon frequencies are very sen- 
sitive to the interlayer distance in graphite [§]. 

Graphite, with its well separated and polarizable 
graphene sheets, is a very good candidate for the ex- 
ploration of LFE-induced phenomena. EELS and IXS 
measurements (see e.g. [H 0-111) have determined the 
energy-loss spectra and plasmon dispersion for a wide 



range of momentum transfer q, and the continuous 
change of the spectrum with change in direction of q has 



been studied [ll|. Calculations of the energy- loss func- 
tion that are based on the homogeneous electron gas (see 



e. g. |12l-ll4|') or the tight-binding model (see e. g. [15l.ll6l|) 
have been used extensively to study in-plane properties 
of graphite. Ab-initio calculations based on Density- 
Functional Theory (DFT) in its time-dependent exten- 
sion (TDDFT) 17], either in the adiabatic local density 
approximation (TDLDA) or even in the Random Phase 
Approximation (RPA), reproduce experimental plasmon 
spectra with very good precision [i,[ll|. These calcula- 
tions have been successfully used to describe the angular 
dependence of energy-loss spectra for relatively moderate 
momentum transfer q. The range of larger q, instead, is 
much less explored. 

In the present work we demonstrate that this range 
offers access to striking phenomena. In particular, our 
ab-initio calculations and IXS measurements reveal an 
anomaly in the angular dependence of the dynamic struc- 
ture factor: For a momentum transfer close to certain 
reciprocal lattice vectors, we observe drastic changes in 
the spectra upon small variations in q. This discontinu- 
ous behavior should have important implications for the 
interpretation of measurements close to Bragg reflections 
in any strongly inhomogeneous system. 

We performed first principle RPA calculations of the 
momentum resolved and frequency dependent dynamic 
structure factor S{q,uj), which is directly related to the 
energy-loss function. The electronic ground state of 
graphite (Bernal stacking) was calculated in DFT-LDA 
(local density approximation) with Abinit using a 
plane-wave basis set [2^ and norm-conserving pseudopo- 
tentials The Kohn-Sham bandstructure was then 

used to compute the independent particle polarizability 
X° with the DP-code [2^. As a next step we calculate 
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the longitudinal dielectric function e and its inverse e~^, 
that links the total to the external potential in linear re- 
sponse, (ptot — e^^<y5oxt- In R-PA t — \ — ux°, where v 
is the Coulomb interaction. In a periodic system e is a 
matrix in reciprocal lattice vectors (G,G'), and a func- 
tion of the reduced component of momentum transfer 
inside the first Brillouin zone and of frequency w. For a 
given momentum transfer q = -I- Go (where Go is a 
reciprocal lattice vector) the dynamic structure factor is 
then (no denotes the average electron density) 



q,=0 



q^=1/8 (-0.37 1/A) 



^(q,^) 



only diagonal elements of are needed. However, as 
e is not diagonal for an inhomogeneous material, its 
inversion will mix matrix elements. One can under- 
stand the physics of this mathematical fact by expanding 
e-' = {l-vX°)-'-- 

+ X]^Go^GoG(qr,w)uGXGGo(qr,w) + ... 
G 

The first order term gives the response of the independent 
particles to the external potential. The second term is the 
response to the Hartree potential that is induced by the 
first order response. This self-consistent process is then 
continued to infinite order. As xo is a matrix, an external 
potential with momentum q^ -I- Go can induce spatial 
charge fluctuations, whose momentum q^ + G differs by 
any reciprocal lattice vector, and to which the system 
will also respond; these are the LFE. 

Most often LFE are studied where the external per- 
turbation is of very long wavelength (q — >• in the opti- 
cal case), and results are modified by contributions with 
shorter wavelength. However, when q = q^-f Go is bigger 
than a reciprocal lattice vector one can also have contri- 
butions from larger wavelength, i.e. smaller |G| < |Go|. 
For example, LFE can couple plasmons from the first 
Brillouin zone q^ with excitations at large momentum 
transfer q = q^ -I- Go, leading to a Fano resonance in 
silicon [20] and lithium 2l| or a periodic plasmon disper- 
sion for momentum transfer perpendicular to the planes 
in magnesium diboride MgB2 [2|. We will show in the 
following that this situation is particularly striking for 
graphite, where it gives rise to an unexpected disconti- 
nuity of the dynamic structure factor. 

Figure 1(a) shows the calculated dynamic structure 
factor S{q,uj) of graphite for a series of q = (0,0,53) 
that are chosen exactly perpendicular to the planes [Note 
that (91,92,93) = gibi -I- (72b2 -I- gabs where the b„ are 
primitive reciprocal lattice vectors; see Fig. 3(b)]. The 
structures that can be seen in this low-energy range are 
only due to the very weak interaction between graphene 
sheets; the intensity of the spectrum is therefore quite low 
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FIG. 1: (color online) Calculated dynamic structure factor 
^(q, oj) of graphite for different momentum transfers q = 
(gi, 0, ga) (a) exactly perpendicular to the planes {qi = 0) and 
(b) with a small in plane component (gi = i). Dashed lines 
indicate spectra for integer ga. All spectra have been shifted 
equidistantly in y-direction for better visibility. Dotted lines 
are guides to the eye. 



and the dispersion small, contrary to the case of MgB2 
[3] where the coupling between layers is much stronger. 

Instead, Fig. 1(b) shows calculations for momentum 
transfers q ~ (|jO, 93) where the in-plane component 
91 = i is kept fixed, and the perpendicular component 
93 is varied. Although for large 93 the direction of q de- 
viates only slightly from the c-axis, we observe a striking 
difference compared to the on-axis results: First, the in- 
tensity is significantly increased for most of the spectra 
- when 93 is zero, we simply see the in-plane 7r-plasmon. 
Second, there is a strong dispersion in the peak positions; 
the main peak shifts between 7.6 eV at 93 = and 6 eV 
at 93 = 1. The most striking observation, however, is the 
disappearance of the peak at 93 = 2: for that value, the 
dynamic structure factor abruptly becomes completely 
flat below 8 eV. There is hence a significant change of 
the spectra for a small change of q. Even more impor- 
tant, the effect cannot be explained by the finite size of 
the in-plane component 91 : in the following we will show 
that these results remain valid for arbitrarily small 91 . 

The bottom panel of Fig. 2(a) shows S{q, uj) calculated 
with and without LFE (solid and dashed lines, respec- 
tively) for momentum transfers q = that are vanish- 
ingly small (|j7| — 5-10^^ A ) and differ only in the 
angle 6 to the c-axis. For in-plane momentum trans- 
fers {0 — 90°), we find a pronounced 7r-plasmon peak 
at 7 eV, whereas in perpendicular direction {9 — 0°) a 
peak, although of much lower intensity, is found at 4 eV. 
For intermediate directions of rj one finds a continuous 
behavior when the angle 9 is changed, as one would ex- 
pect. LFE contribute only marginally, i.e. the anisotropy 



3 



2 q^0,0,2) + Ti 




p q=(0,0,0) + T1 

x5 





q = + Go 



(a) 



2 4 6 8 10 

energy loss (eV) (b) 



FIG. 2: (color online) (a) Dynamic structure factor S{ri + 
Go, oj) for small deviations |t7| = 5 ■ 10~^ A ^ from reciprocal 
lattice vectors Go calculated with (solid) or without (dashed) 
LFE. (Bottom) For Go = (0,0,0), S{-n,uj) depends on the 
direction of rj. (Middle) For Go = (0,0, 1), arbitrary r] lead 
to the same spectrum. (Top) For Go = (0, 0, 2), S{r] + Go, oj) 
again changes with rj, but only when LFE are included, (b) 
Definition of the angles 8 and 6' . 



simply stems from the anisotropy of the head element 
eoo{ri,Lu) of the dielectric matrix. These results are con- 
sistent with earlier calculations on graphite 

Moving on to large momentum transfers, the next 
higher panel displays the results for q = (0, 0, 1) + rj. A 
change in rj corresponds now to an infinitesimal change 
in q [see Fig. 2(b)], and results both with and without 
LFE do not vary with rj; therefore, only one result is 
shown. It should be noted however that LFE start to 
become significant, because of the increased momentum 
transfer. 

This picture changes completely when we reach mo- 
mentum transfers q = (0, 0,2) +r] near the second recip- 
rocal lattice vector, shown in the top panel of Fig. 2(a): 
whereas the spectrum without LFE is completely flat and 
remains stable while t] is changed, we flnd that LFE 
lead to drastic modifications of the spectra for infinitesi- 
mal variations of the total momentum transfer q. A di- 
rect comparison between spectra for momentum transfers 
q = (0, 0, 2) -|- ?7 and the corresponding reduced momen- 
tum transfer rj shows that they are very similar besides 
a scaling factor, whenever r] is not exactly in-plane (top 
vs. bottom panel); in other words, LFE lead to the reap- 
pearance of spectra of lower Brillouin zones 0, HO] , and 
the dynamic structure factor S near the reciprocal lat- 
tice vector (0, 0, 2) is determined by the direction of the 
reduced momentum transfer q^ (here q^ — r)) and not by 
the direction of q itself. 

Graphite is a very convenient case for further analysis. 
Indeed, in the case of small q^ and Go = (0,0,2), we 



find a particularly strong coupling between q^ and q^ -I- 
Go as the wing element eoGo(qr,w) dominates by far all 
other off-diagonal elements eGG'(qr,w). Consequently, 
neglecting the coupling to all other modes, the inversion 
of e reduces to the inversion of a simple 2x2 dielectric 
matrix [23] and one obtains (here, the w dependence is 
omitted): 



,(qr) 



1 



eGoGo(qr) 



£Goo(qr)£OGo(qr) 
[eGoGo(qr)]2 



10 i^r 



This result is similar to the two-plasmon-band approx- 
imation (22j . The first term corresponds to the result 
without LFE. The second term leads to the reappear- 
ance of the spectrum eQQ(qr,CLi). The latter is however 
weighted by the off-diagonal elements. Whenever the sec- 
ond term is strong, the spectrum for large q = q^ -I- Gq 
will also depend on the one for the reduced component 
qr, and hence on its anisotropy. 

Still, it remains to be understood why this coupling to 
q — >■ does not show up neither for in-plane deviations 
q = {i], 0, 2), nor around Go — (0, 0, 1). To this end, we 
can make use of general properties of the dielectric func- 
tion in semiconductors. In the limit of high frequencies, 
the coupling element eoGo (qr, w) can be approximated as 



I , 47r q^ 
eoGolqr,^^) = ^ — 



Gn 



ri(-Go), 



where n(— Gq) denotes the Fourier coefficient of the elec- 
tron density. A similar expression has been found by 
Sturm and Oliveira [2^ in the framework of a quasi- 
free electron gas. First, we see from this equation 
that the coupling between q^ -I- Go and q^ is propor- 
tional to the cosine of the angle B' between the two 
corresponding directions [see Fig. 2(b)]. The prefactor 
/ = EGooeoGo/eGoGo °^ ^'^^ ^' enforces the anisotropy of 
the spectra: in particular, for a small in-plane q^ one has 
9 = 90° K, 6' , which explains the absence of strong LFE in 
the spectrum for 6 — 90° [Fig. 2(a), top panel]. Second, 
the coupling vanishes whenever the Fourier component 
n{—Go) of the electron density becomes zero. As this 
coefficient is proportional to the crystal structure factor, 
a wing element vanishes whenever the Bragg reflection of 
the corresponding reciprocal lattice vector is forbidden; 
for graphite in Bernal stacking, this is the case for all 
Go ~ (0, 0, 2m + 1) where m is an integer. LFE around 
(0, 0, 1) stem hence only from a mixing with other non- 
vanishing G 7^ beyond the 2x2 model. They do not 
introduce any signiflcant dependence on the direction of 
qr as G -I- qr w G. Instead, for Go — (0, 0, 2m) the two 
graphite planes A and B in the unit cell contribute in a 
constructive way (analogous to the constructive interfer- 
ence in the case of the Bragg reflection), which leads to 
the strong effect. 

With these explanations in mind, let us come back to 
the results shown in Fig. 1(b). Since q,. is still reasonably 
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FIG. 3: (color online) (a) Comparison of the structure factor 
^(q, oj) measured by IXS-experiments (dots) and calculated 
in RPA (lines; dashed for integer 53) for q — (i,0, ga). The 
elastic tail has been removed from the raw experimental data 
and a uniform scaling has been applied, (b) Measured q- 
points in reciprocal space. 



small, the above arguments hold. In particular, we can 
explain the drastic spectral changes near (0, 0, 2) by the 
fact, that (i) the spectra from the first Brillouin zone 
reappears, which strongly depends on the direction of 
(angle 9) due to the anisotropy of graphite, and (ii) the 
coupling and hence the strength of the recurring spectra 
is proportional to |qr -ql^ (angle 9'). 

One may wonder, whether such a spectacular effect can 
actually be observed, or whether it is masked either by 
exchange-correlation effects that are neglected in RPA, or 
by experimental conditions like a strong elastic tail near 
the Bragg reflection, that prohibits a direct observation 
of the discontinuity at (0,0,2). Therefore, we have con- 
ducted measurements of inelastic x-ray scattering for mo- 
mentum transfers q — (i,0, ^) (for selected numbers n 
between 6 and 16) covering the range between and above 
the Bragg reflections (001) and (002). The in-plane com- 
ponent Igibil = 0.37 A is still small enough to create 
the desired effect [see Fig. 1(b)], but large enough to avoid 
the Bragg reflection. The measurements were carried out 
at the Taiwan inelastic scattering beamlinc in SPring-8 
(BL12XU). The graphite sample was a plate having a 
surface of 2 x 3 mm^ and a thickness of 100 \im 2J|. The 
x-ray Laue photograph showed very nice spots, indicat- 
ing that the sample was not Highly Oriented Pyrolitic 
Graphite but a single crystal. The energy resolution was 
140 meV. A Si 400 four-bounce monochromator and a 
Si 444 spherical crystal analyzer were used. The momen- 
tum resolutions were approximately 0.15 A along the 
horizontal and the vertical axes. In order to subtract the 
tails of the elastic lines for the spectra, glass was also 



measured as a reference. Fig. 3(a) shows the measured 
curves (dots) , together with the corresponding calculated 
results (lines). The agreement is very good; in particu- 
lar, the predicted peak shift is clearly seen in the mea- 
surements, as well as the abrupt change from a peaked 
spectrum for q = (i,0, -^) to a completely flat one at 
q = (|,0, 2), and the difference between = ^ and 
due to the different angles 9' [see Fig. 2(b)]. Hence, 
our measurements give unambiguous support to the pre- 
sented theoretical prediction. 

In conclusion, our RPA calculations and IXS measure- 
ments revealed and explained a striking discontinuity in 
the dynamic structure factor S'(q, lo) of graphite at mo- 
mentum transfers q^ -f Gq close to Bragg allowed recip- 
rocal lattice vectors Gq: infinitesimal changes in the mo- 
mentum transfer induce strong changes in the resulting 
spectra. No discontinuity is observed when the crystal 
structure factor vanishes. It is hence a consequence of the 
Bernal stacking of the graphene layers that no changes 
occur at Go — (0,0, 1). Such a behavior of 5(q, w) has 
important consequences for theory and measurements: 
First, anisotropic excitations from the first Brillouin zone 
might appear at large momentum transfers q leading to 
an anomalous angular dependency of the spectra. Sec- 
ond, these recurring excitations belong to a direction q^ 
different from q. Third, from the experimental point of 
view, whenever measurements of the dynamic structure 
factor are performed close to an allowed Bragg reflec- 
tion, the resulting spectra can be extremely sensitive to 
the chosen momentum transfer. We expect similar ob- 
servations in other anisotropic crystals that show strong 
crystal local field effects, especially for layered or quasi 
ID structures. 
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